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Abstract 

We use the relation between 2d Yang-Mills and Brownian motion to show that 2d Yang-Mills 
on the cylinder is related to Chern-Simons theory in a class of lens spaces. Alternatively, this 
can be regarded as 2dYM computing certain correlators in conformal field theory. We find that 
the partition function of 2dYM reduces to an operator of the type U = ST P S in Chern-Simons 
theory for specific values of the YM coupling but finite k and N. U is the operator from which 
one obtains the partition function of Chern-Simons on S 3 /Z p , as well as expectation values 
of Wilson loops. The correspondence involves the imaginary part of the Yang-Mills coupling 
being a rational number and can be seen as a generalization of the relation between Chern- 
Simons/WZW theories and topological 2dYM of Witten, and Blau and Thompson. The present 
reformulation makes a number of properties of 2dYM on the cylinder completely explicit. In 
particular, we show that the modular transformation properties of the partition function are 
intimately connected with those of affine characters. 
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1 Introduction and summary of results 



Two-dimensional Yang-Mills is a well-studied theory [1, 2, 3], especially in its large N limit 
where it can be reformulated as a string theory [2] . Being a nice playground to study quantum 
gauge theory, both perturbatively and non-perturbatively, as well as being able to explicitly 
manifest itself as a string theory, this theory very explicitly realizes the idea of 't Hooft [4] that 
gauge theories at large N can be reformulated as string theories, and has been an object of 
extensive study for the last twenty years or so. To the extent that one might think there are no 
interesting questions for Yang-Mills left to answer on a Riemann surface. 

Yet there are several indications that much physics is still to be extracted from bosonic 
2dYM. Let us mention just a few of them. In [5], it was shown that 2d Yang- Mills on the torus 
is related to a string theory at large N, but of a different kind from the one found in [2, 6]. It is 
not defined on a two-dimensional target space, but rather as the topological A-model on a six- 
dimensional non-compact Calabi-Yau. The correspondence was argued both on the topological 
vertex side and on the 4d SYM theory side, where the twisted 4d SYM reduces to 2d bosonic 
YM. Related ideas had been discussed previously in [7]. Vafa has further conjectured [5] that 
2dYM may have some relevance in the description of 4d black holes of type II superstrings. 

Nekrasov [8] has shown that random partitions can be used to compute instanton contribu- 
tions to the prepotential of = 2 SYM theory. In [9], it was shown that these partitions are 
related to, among other things, Wilson's plaquette model of QCD2 at finite N. In [10], random 
partitions were shown to be related to large N 2dYM on the sphere. 

In 2dYM, the partition function of the sphere and the torus can be obtained from that of the 
cylinder by sending the external states to zero and by tracing over them, respectively. Therefore, 
the above suggests that a reexamination of 2dYM on the cylinder might be desirable. 

Finally, in [9] a relation was found between Brownian motion in a Weyl chamber and Chern- 
Simons theory, and between Brownian motion in a Weyl alcove and 2d Yang-Mills. This sug- 
gested some type of relation between Chern-Simons and 2d Yang-Mills. Such relations had 
already been pointed out in [11, 12, 13]. 

In this paper we reformulate 2d Yang- Mills on the cylinder along the lines of [9]. We use 
the techniques of earlier studies of the heat equation in the mathematical literature [14, 15]. 
This makes the modular properties of the theory completely explicit. In [16, 17], these modular 
properties were studied for the torus. The chiral theory was found to transform with weight 6g— 6 
in the 't Hooft expansion, but the non-chiral theory did not seem to have definite transformation 
properties. On the cylinder, we can reformulate the partition function of 2d Yang-Mills in terms 
of generalized theta-functions and work out its modular properties explicitly. 

The new formulation of 2d Yang-Mills also sheds some light on the relation between 2dYM 
and Chern-Simons. We find that for special imaginary values of the coupling (2ir times a rational 
number) and particular external states, corresponding to integrable representations, the partition 
function equals the gluing operator U = ST P S in Chern-Simons theory, from which the partition 
function on S 3 /Z p and expectation values of up to two Wilson loops can be obtained. Notice 
that the partition function of 2d Yang-Mills diverges at these particular values of the coupling. 
This divergence is analogous to those found in [12, 18] at zero coupling, and reflects the fact 
that 2dYM is generically non-analytic in the coupling. 

In another line of developments, Witten [11, 12] and Blau and Thompson [18] showed that 
topological 2dYM/BF theory is related to Chern-Simons theory on certain spaces, and in partic- 
ular computes the dimension of the space of conformal blocks in the WZW model. Topological 
2dYM is simply the g Y M — > limit of 2dYM, where the partition function reduces to the volume 
of the moduli space of flat connections on a Riemann surface. It is equivalent to BF theory with 



2 



vanishing potential. This can be compared to CFT only in the limit k — ► oo. To go beyond 
the k — > oo limit, Blau and Thompson [18] showed that one can use the same topological BF 
theory but with a compactified scalar. The limit k — > oo then corresponds to the limit where the 
scalar decompactifies. This theory however did not appear to have a straightforward Yang-Mills 
interpretation for finite k. 

Our results could be interpreted as a 2d Yang-Mills description of the finite k case, giving 
some sort of topological sector of 2dYM for special values of the coupling. Our theory can be 
reformulated as an ordinary BF theory (i.e., with non-compact scalar) with a quadratic potential 
and complex coupling. 

The organization of the paper is as follows. In section 2 we analyze the structure of 2dYM 
on the cylinder by reformulating it in terms of Brownian motion and obtain the expressions that 
are used in the remainder of the paper. A reader who is mainly interested in the relation to 
Chern-Simons/WZW model might want to have only a quick glance at this section. In section 
3 we show that complexifying the Yang-Mills coupling the partition function of 2dYM is equal 
to the (appropriately normalized) operator U = ST q S. In section 4 we discuss how our results 
relate to those by Witten, and Blau and Thompson where they compute dimensions of moduli 
spaces of flat connections by using Chern-Simons theory. We end with some conclusions and 
open problems. 

2 Two-dimensional Yang-Mills on the cylinder 

In [9] a relation was found between Brownian motion/discrete random walk models of N non- 
intersecting movers on a line, and certain low-dimensional gauge theories. The gauge theories 
involved were 3d Chern-Simons theory, 2d Yang-Mills, and the plaquette model of QCD2. As it 
turned out, a fruitful way to think about these motions is in terms of a single particle performing 
Brownian motion in the Weyl chamber of some (finite or affine) Lie algebra. In this section we 
will briefly review the basic correspondence and work out in detail the correspondence between 
Brownian motion and 2d Yang-Mills, leaving the Brownian motion description of Chern-Simons 
theory for a future publication. 

Let us first recall some basic facts about 2d Yang-Mills on the cylinder that we will use later 
on. The partition function can be shown to be [1] 



The Yang-Mills coupling always appears in combination with the area of the cylinder, t = fl^M^- 
C 2 {\) is the eigenvalue of second Casimir in the representation A (for notation and conventions, 
see appendix A). The sum runs over all finite dimensional, irreducible representations, labeled 
by their highest weight. 

This partition function is the kernel for the heat equation on the group manifold of the 
corresponding group G: 




(1) 



AeP+ 



-KiZ 2 dYM{g,g') = -77 C 2 (g)Z 2dYM (g,g') ■ 



(2) 



Given a boundary condition fo(g) at t = 0, the heat equation has a unique solution: 




(3) 
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We will later on see that it is important that the coupling constant t = <7ym^ * s complex, and in 
particular not purely imaginary. The reason is that the partition function is not analytic when 
the real part of the coupling goes to zero [12]. Thus, one should be careful when continuing 
solutions of the heat equation to solutions of the Schrodinger equation because in general the 
continuation is not analytic. We will show this explicitly. 

The large iV limit of 2dYM on the cylinder has been studied in [19]. We should here also 
mention the interesting related work [20]. 

2.1 Brownian motion and low-dimensional gauge theories: the basic corre- 
spondence 

In [9] we considered Brownian motion of a number r non-intersecting movers on a line. Alter- 
natively, this can be regarded as motion of a single particle in the Weyl chamber of a Lie group, 
whose dimension equals the rank of the gauge group, r. We will first review the prototypical 
case of Chern-Simons on S 3 with gauge group U(N), in which case r = N. In what follows we 
briefly summarize this relation [9]. Details will be presented elsewhere. 

Even though for some purposes the description in terms of N particles is more useful than one 
where one considers a single mover, in this paper we will describe things from the perspective of 
a single mover performing Brownian motion in the Weyl chamber of some Lie algebra. Consider 
first U(N). Its Weyl chamber is the dual Cartan subalgebra modded out by the Weyl group, in 
this case the permutation group Sn. We can arrange coordinates so that the Weyl chamber is 
given by a vector A = (Ai, . . . , Ajv), where Ai > A2 > . . . A^. The probability for this particle to 
travel from p to A in time t is given by [21] 

Pt ' N{Xl / " ) = (27rtV /2 e ~^^ det \ eXiN/t U<i<j<N • (4) 

The vectors A and p label irreducible representations of U(N) when they lie in P + , i.e. when 
their coefficients are positive integers in a basis of fundamental weights, which we will assume 
in what follows. Notice, however, that Brownian motion is a continuous process and so the 
intermediate steps do not lie on this lattice, but only inside the Weyl chamber. Now it was 
shown in [9] that this is the partition function of Chern-Simons on S 3 with gauge group U(N) 
when the particle travels around a loop starting and ending at A = p = p. p is here the Weyl 
vector, denoting the trivial representation, and plays the role of our origin of coordinates on 
the Lie algebra. By taking different boundary conditions one was also able to get expectation 
values of Wilson lines, and precise agreement was found. As noted in [9], the above probabilities 
together with the Boltzmann factor provide a representation of SL(2, Z) on affine characters and 
so for boundary conditions corresponding to integrable weights they are in 1-1 correspondence 
with the modular matrices of the WZW model. More generally, they are characters of the 
corresponding gauge group. 

An alternative expression for the probability was given: 

which is obtained by the method of images [21]. 

2.2 2d Yang-Mills as a Brownian motion 

In this section we consider a particle moving in the fundamental Weyl chamber of the affine Lie 
algebra of the corresponding group. The structure of this Brownian motion is quite interesting, 
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and it turns out to be related to 2d Yang-Mills on the cylinder. 

In appendix A we recall some definitions and notation that we use in what follows. 

We will be concerned with Brownian motion in the fundamental Weyl alcove of the algebra. 
This is the set of weights, moded out by the affine Weyl group W. We recall that the affine 
Weyl group also includes translations in the coroot lattice Q y : 

W = W k T , (6) 

where W is the Weyl group. Thus we can characterize the fundamental Weyl alcove as P/IQ W , 
where I = k + g - here, k is the level and g is the dual Coxeter number -, moded out by 
the Weyl group. Notice that the Weyl alcove is actually finite, so we can think of the non- 
intersecting particles as moving on the compact dual Cartan torus. Therefore, in order to 
obtain the Brownian motion density in the Weyl alcove we sum over all the images under this 
translation [9]: 

MA.M) = jA^ E E e(w)e-^-^ 2 (7) 

where 1 |A| 2 = Yli=i A?- F° r U(N), r = g = N and I = k + N, but in the following we will not 
assume this. 

By construction, <7t, r (A, /x) still satisfies the heat equation 

d 1 

— qt, r (\n) = j A 9t,r(A,/x) (8) 

where A is the flat space Laplacian with respect to the boundary conditions A or fi. 

There is a different way to rewrite the above expression which is in some respects more 
natural. It is not hard to see that the following identity holds: 

E eMe^ l7+A ^ 12 = e'^+W 2 ) £ e (w)e> X '^ (9) 

where we defined affine vectors 2 A = A + ICoq, jl = \i + Iujq. Notice that these vectors are in the 
fundamental alcove of the algebra, as they should. In fact, they are of the form A + p, where A 
has level I. Thus we get 



1 

r/2 



<//./(A./0 =^72 E • (10) 



We could have guessed the above expression from the analogy with the finite case. In fact, 
discretized versions of such expressions for the case of random walks are mentioned in [22, 23] 
and analyzed in [24] in the context of a generalized ballot problem. For Brownian motion, [25] 
gives a related expression for motion of iV non-interacting particles on a circle. The above gives 
a derivation based on the compactification of the infinite Weyl chamber. 

We now turn to the relation of the Brownian motion density to 2d Yang-Mills theory on 
the cylinder. Our aim is to show how the partition function of 2dYM, (1), is related to the 

1 Despite the notation, the absolute value bars do not imply taking complex conjugates, but are simply a 
Cartesian product in a suitable basis. 

2 We can include the quadratic factors in the definition of the affine weights by adding non-zero n\ pieces. 
However, we prefer to keep these factors explicit. See the next section. 
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Brownian motion density, (9). To that end we make use of Frenkel's results [14]. We rewrite 
theorem (4.3.4) of [14] as follows: 

where D p {\) is the denominator of the character: 

D p (X) = Yl 2sinh(a,A)/2 . (12) 

P is the weight lattice, and \P/Q V \ is the volume of P/Q y , i.e. the number of points of P that 
are contained within a unit cell of the coroot lattice Q v . For the external states on the cylinder 
we filled in group elements g' = e A , g = e M , and we defined 

A' = ^A 



I 

2m 



» = — H 

The proof, which we will not repeat here, goes in two main steps. First one uses Weyl's character 
formula on the left-hand side and reduces the final expression to an expression involving only a 
sum over weights and a single sum over elements of the Weyl group. Then one uses Poisson's 
resummation formula to replace the sum over the weights by a sum over coroots. 
Using this, we find: 

< 14 > 

This is the expression in [9] . Notice that the Poisson resummation effectively transforms the 
time-dependence as t — > — 1/t. 

Notice that the relation between 2d Yang-Mills and Brownian motion involves factors of 
Dp (A) that depend on the representation. These come from the normalization of the characters 
X\(g)- The partition function has been normalized so that it has the expected limit on the disc, 
i.e. when g — > id. Of course, such normalization factors are independent oft. This normalization 
ensures that both sides satisfy the heat equation and will be explained in section 4. On the other 
hand we see that Chern-Simons and Brownian motion share the same normalization. 

2.3 The structure of 2d Yang-Mills on the cylinder 

In [9] we claimed that the modular properties of 2d Yang-Mills on the cylinder are simple to 
analyze. This follows from the representation (14) and the analysis is absolutely standard, see 
e.g. [26, 27]. Since the normalizing prefactor has no ^-dependence, we will analyze in detail the 
expression 

= — l -j- 2 £ £ eMe-^-^l 2 . (15) 
This can be expressed in terms of the following theta-function: 

e(J)(£. r . s ) =e -2irils £ e rfr|a v +A/«-C/r| 2 -rf|C| 2 /r ( lg ) 
o v GQ v 



6 



where I is any integer (in our case it will be the level). Taking 

1 2-Ki 

t = T 

r 

. = -f. <») 

we have to evaluate the 0-function at the point (£; t; s) = — ^ — ^r) ■ Thus we get the 

following identity 

qt,r(\,n) = {2iTt y /2 E e ( w ) Q w(\)(C,r;s) . (18) 

^ ' w€W 

S'-transformations are generated by (£; r; s) — > — ^; s + . In this case, the S'-matrix is 
represented by 

S (0 = ^Zll |p / QV r l/2 e -2 7ri(M ,A)// (19) 

and 

e ?(;'4' + ! 5r) = S S«e!?«;r;»). (20) 
V J veP/lQv 

The sum 

<MC;r; S )= E £ W4(C; t ;') ( 21 ) 

transforms under ^-transformations as: 



^(^~;* + ^) = (-^) r/2 |P/QT 1/2 E e H E e-^^eW^r;-) (22) 



and under T-transformations: 

A (C;r + l;5)=e i ' r W a /^ A ( C;r;5 ) . ( 2 :;) 



It is now straightforward to see how qt, r (\ A 4 ) transforms. First of all, let us add the additional 
parameter s in our notation. This simply keeps track of the phases generated: 

q t>r (\,H;s) = e- 2mls q t>r (\^) • (24) 
A modular transformation on (C;r) — > (C'jt') redefines fi,t as follows: 



1 [2m\ 



2 

t 



H' = /- . (25) 



We get: 

r/2 



M^'iO = (j)' E £ H E ^L e ?(C;r; S ), (26) 
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where S is defined as before. This can also be written as 



= i-^\+{2Ktyl 2 Y, S Xv q t , r {v,^s) (27) 
uePl + 

where S is now the usual, unitary representation of the S'-matrix on affine characters. Notice 
that in our case we start with the particular value of s (17) and end up with s' = 0. Now because 
of (14), this gives us the transformation rules of Z 2 dYM(9, g'\ t). We see that qt,r transforms in 
a unitary representation of 5L(2,Z), up to a weight of r/2. 

We get a unitary representation of the modular transformations if we appropriately normalize 
the partition function. Indeed, it should by now be clear that 

= e ~^ (A) *xW) (28) 

Qt,r(P, M) 

where the affine character is 

A and fi are defined as in (9). In our case the modular anomaly term is absent because 3 n\ = 
ria = 0. The above formula is completely analogous to Brownian motion in the Weyl chamber 
of the classical group, for example U(N), where we found that the normalized probability was 
the character of U(N). The modular transformation properties of the affine character are now 
the usual ones: 



x -(C;t + 1;s) = £ T x . XA (C;r; S ) 
Aep£ 

^( C r'-h S+l €) = E S AA^(C;ra), (30) 

AeP£ 

where S and T are the usual representations of the modular matrices on affine characters, given 
in appendix A. 



3 Two-dimensional Yang-Mills on the cylinder and C hern- Simons 

on S 3 /Z q 

Some experience with the modular matrices of the WZW model and U(N) characters imme- 
diately tells that, up to a normalization, an object such as (1) is very similar to a product of 
modular matrices ST q S for some power q related to the YM coupling. There are, however, 
important differences. The first one is that the Hilbert spaces are different. In 2d Yang-Mills, 
we sum over all finite dimensional representations, and also the external states at the two ends 
of the cylinder can be in arbitrary representations. The WZW model, on the other hand, is 
based on an affine Lie algebra where the Hilbert space only contains integrable representations 

3 It will give a non- vanishing contribution, however, if we decide to include the Casimir in the definition of the 
affine weights, as discussed earlier. 
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satisfying (A, 9) < I, where 8 is the highest root. This Hilbert space is finite. The second differ- 
ence is in the coupling constant. 2d Yang-Mills really satisfies a heat equation, whereas ST q S, if 
q or the level were continuous parameters, would satisfy the Schrodinger equation. Nevertheless 
one would expect that some kind of analytic continuation should be possible, in the same way 
that we can describe topological strings by identifying the string coupling as g s = 

Of course, we expect that in the large k limit, and up to possible divergences, 2dYM and 
Chern-Simons agree, as in that case the Hilbert spaces are the same. This is indeed the result 
of [11, 28, 18]. 

The purpose of this section is to make the above observations on the finite k case precise. To 
this end, it will prove far more useful to start from the expression (14) than to work directly with 
(1). Basically, in (14) the two characters have already been expanded into a single expression of 
the type found in the Chern-Simons literature. 

It turns out that there are special values of the coupling g YM at which the partition function 
(1) collapses to the gluing operator ST g S in the way expected, i.e., up to normalization. These 
values of the coupling are imaginary, in particular they are 2iri times a rational number, and 
precisely at these values the sum over all representations reproduces the sum over the integral 
representations in Chern-Simons. This uses a certain periodicity of the partition function for 
these values of the coupling. Since the partition function is periodic on an infinite lattice, there 
is an overall diverging factor multiplying the partition function that we will explain. 



3.1 Chern-Simons on S 3 /Z q 

In [28] , a surgery prescription was constructed for obtaining partition functions and expectation 
values of Wilson loops on non-trivial manifolds. One glues together a three-manifold M out 
of two solid tori by appropriate identification of the A- and -B-periods of the two tori after 
diffeomorphisms. These diffeomorphisms are represented by SL(2, Z) transformations. The 
partition function of Chern-Simons theory on this manifold is then given by the representation 
of the corresponding gluing operator on the Hilbert space of the WZW model. 

Chern-Simons theory in lens spaces was introduced in [28] and studied in [29, 30, 31, 32]. 
The corresponding matrix model for the case S 3 fE q was analyzed in [32, 33]. It was shown in 
[32] that modding out the S 3 by 7L q corresponds to acting on the cycles of the torus before gluing 
by an SX(2,Z) matrix U = ST q S. If we represent a generic SX(2,Z) matrix by 

U(P,q) = ( P r 

\ q s 

and in the usual representation 

s- (;-')■ <»> 

then it follows that we need 4 U^ 1 '^ = ST Q S with p = s = — 1, r = 0. Flipping the signs of p 
and s would simply correspond to flipping the sign of g YM and of the external states. We will 
not consider this in this paper. 

4 In this, our conventions differ from those in [32], where all signs in U were positive. Of course this is a trivial 
redefinition, as one can explicitly show. In fact one matrix is minus the inverse of the other. 
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The representation of the generic operator U^'^ in this Hilbert space is given in [29, 30, 31]. 
In the case of interest to us we have: 



(_!,,) _ [*sgn(g)]l A +le g^g^j ^ lh+x+Wfl{ 2 

^ ~ (l\q\)r/i lP/QVli/2 Z. Z. e W • ^ 

Here, ^(i/fr' 9 )) is the Rademacher function, related to a choice of framing [34, 12, 29]. In this 
paper we will not keep track of constant phase factors, although this would be an important 
check. 

As shown in [30], taking matrix elements of U^ 1 '^ in the trivial representation gives the 
partition function of Chern-Simons on S 3 /Z q . Taking it between non-trivial representations A 
and \i means that we glue two tori with colored cycles A and fx, which results in the expectation 
value of a link W\ p embedded in S 3 /Z q . 

As a trivial check, the geometric picture tells us that the case q = 1 should give us back the 
three-sphere. This is indeed the case, as we can explicitly see from (33), which (up to framing) 
reduces to T~ 1 ST~ 1 . This is nothing else than STS, as it should. 

We can already see the very similar structure with (7). In fact, under the identification 
1/t = 2iri/lq, up to numerical factors independent of I and q the only difference seems to be 
the summation range. One would think that in the limit q — > oo both expressions precisely 
agree. Things are nevertheless a bit more subtle. The expression (7) does not make sense when 
1/t = 2m/lq. Due to the periodicity of the lattice, an infinite number of images will contribute 
exactly the same weight. That is why, whereas the real expressions are finite, continuation to 
complex values of the coupling is subtle in this case. In fact we will exploit this periodicity 
to show that, up to this divergence, both expressions do agree without the need to take q to 
infinity. The divergence is well-known in 2dYM. The partition function of 2dYM is in general 
non-analytic as g YM — ► (depending on the Riemann surface and gauge group) [12]. We seem to 
find here that it stays non-analytic for these special values of the imaginary part of the coupling, 
whenever Re<? Y M = 0. Although from the point of view of 2dYM it is not always natural to 
remove these divergences [11, 12], they can be consistently absorbed in the definition of the 
zero-point energy, as we will discuss later. 

Prom the Brownian motion point of view, we might say that complexifying time as above 
effectively amounts to putting periodic boundary conditions inside the Weyl chamber, so that 
we are restricting to the integral weights. All the paths outside this domain then contribute 
with the same weight. 

3.2 Relation between 2d Yang-Mills and Chern-Simons 

Let us now see how we can continue the partition function to imaginary values of the coupling. 
This has to be done carefully, as it diverges on the imaginary axis. We can rewrite the partition 
function as 



V 

1 1Y /2 

27rt) D p (X')D p {-n 



rt) D p (\>)D p (-»>)^^ W)e 



2\ r/2 



\27Tt) D p {\')D p {-n') ^ w { ) 
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where for convenience we included an extra factor of \p\ 2 in the exponent, and as before A' = 
2iri\/l, ix' = 2iri[i/l, and t' = (2iri/l) 2 t. We also defined the function 

z(A) = e~^ 2+2 ^ . (35) 

This is a higher-dimensional generalization of the usual theta function, whose asymptotic prop- 
erties are well-known. It is well-known to converge when i/t is in the upper-half plane, and for 
any A £ C. 

We will be interested in the asymptotic properties of this function. To get the partition 
function from this we just multiply by the relevant Gaussians and sum over images. Since we 
want to compare with Chern-Simons theory, we need to identify 

MM ») 

where e > is a real constant that we will take to zero at the end of the computation. We can 
rewrite the above as 

*(A) = £ e~T (37) 

7e<2 v 

We now split the sum as follows: write 7 = q v + qf3 v where a v 6 Q v /qQ v , /3 V £ Q v , so 

z(\) ^ J2 expH7r(/a v2 /g + /g/3 V2 + 2/a v /3 v + 2a v A/g + 2/3 v A) + 

a v e<3 v /gQ v /3 v GQ v 

-e(/a v2 + 2a v A + g 2 /3 v2 + 2q(3 v (a v + A))] . (38) 
Now the following conditions hold 5 : 

7 2 G 2Z+ 

7A G ZZ . (39) 
We can use this property to reduce the sums. The only terms that remain are: 

Z ( X ) = Yl e" e[9 ^ V2+2 ^ V(QV+A)1 "^ ( ' QV2+2aVA)_e( ' aV2+2QVA) . (40) 

Obviously, the last term proportional to e we can drop, as it has no /? v -dependence and only 
gives a finite contribution that will go to zero when we take e — > 0. The second term in the 
exponential is the one we are interested in; it is the one related to U. The first term is there only 
to regularize the sum and contributes an infinite normalization constant. Clearly, the diverging 
contribution comes from the infinite sum over Q y in the limit e — > 0. The contribution of this 
term is analyzed in the appendix B. Using those results, we finally get: 

^-{^Tm^ £ ■ (41> 



5 Notice that [14] defines the weights with a factor of 2ni, which explains our normalization A' = 2mX/l. The 
rescaling by I was chosen for convenience. 
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Filling this into the summation in (34), we get for the leading divergence 



E £ eWe -A^-P = ( * »^ E E ^)e-^ + ^' . (42) 



The last term is the gluing operator (33), where the irrelevant sign difference in fi comes from the 
fact that at one end of the cylinder we are considering the conjugate representation g -1 = . 
Filling this into the partition function (34), we get 

This is now, up to normalizations, the gluing operator U^' q . As mentioned at the beginning of 
this section, the normalization by factors of D p has to do with the normalization of the character 
in 2d Yang-Mills. This is different from the normalization of S, which does not contain that 
factor in the Weyl denominator formula. If we normalize the gluing operator the same way as 
the characters, that is, such that the S'-matrices are normalized to unit strenght, we get 

«*. = £p§± = (Vkir /2 ^^f;> /Q V |^ E E • (44) 

Thus, we finally obtain: 

/ 2 \ r/2 

u *" = e< *(vJ Z 2dYM(e A ',e^;t') (45) 

where is a phase that we are not considering in this paper 6 . The quantity on the right-hand side 
is finite. The fact that the normalizations now agree precisely matches our naive expectations. 

Coming back to the case of S* 3 , it is now clear that the transformation t — ► — 1/t in (11) 
means that STS = T~ 1 ST~ 1 , as we see from the above. 

The divergence is not new and can now be absorbed as a constant in the definition of the 
zero-point energy. Notice that it depends only on q and e but not on any other quantities. 
Recalling that the action was 

S 2 dYM = -A— I d^v^TrF 2 , (46) 
we can add a counterterm given by 

S ct = cJ^d 2 x^ (47) 

and we see that we need to adjust c = ^ log ^ = Iqr^f-log^. Such counterterms are 
invariant under area-preserving diffeomorphisms and hence respect the symmetries of the theory. 
In our case, the area drops out in (47), and the counterterm is really a pure number, because e 
is. Both g YM and A are dimensionful, but their product is dimensionless: 

9 , 2-iri 2e . . 

9LA = li ~J- (48) 



6 This phase is related to the choice of framing in Chern-Simons and would be interesting to compute. 
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Thus, although we have fixed the product, we still have the freedom to rescale — > <?ym*> 
A — > .A/i, thus preserving the symmetry. 

Counterterms that redefine the zero point energy and thus shift the partition function by a 
constant have been discussed in [11, 12]. There, it was discussed that there is an ambiguity in the 
definition of the partition function related to the addition of such terms 7 , different regularization 
schemes giving different values of c. This divergence is however just a manifestation of the non- 
analytic structure of 2dYM. 

Notice that the divergent prefactor of (e/7r) r / 2 in (45) is precisely the factor of (2irt/l) r / 2 that 
we got in our final expression for our modular transformation (27) if we subtract the imaginary 
part of t. Also, notice that shifts of r — > r + 1 correspond to shifts of t by 2ni/l. This suggests 
that the factor that we found by explicit computation might be related to the prefactor that 
follows from the modular transformation properties, although we have not checked this explicitly. 
In particular, the appropriate g-dependence should be taken into account. It might be possible 
to use the modular transformations in section 2.3 to get the above results directly. 

Let us comment on the normalization of the operator that we used in (44) . We will see 
that this normalization agrees with the results in [11, 12]. Whereas the normalization in 2dYM 
is different from the normalization of U^~ 1,q \ Brownian motion actually automatically comes 
with the same normalization as U^~ 1,q \ This was expected too, as it was found in [9] that, up 
to framing factors, these normalizations generically agree in the case of Brownian motion in the 
Weyl chamber of the finite algebra, S having the interpretation of Brownian motion densities. 
In particular the above implies that one should be able to obtain qt, r {\ A 4 ) from a composition 
of two finite probabilities pt, r with a Boltzmann factor, and summing over intermediate states. 
This should be most clear in the matrix model representation. 

4 2dYM, the WZW model, and the Verlinde formula 

We will now comment on the results in the previous section from a different point of view, 
mentioned in the introduction. In a beautiful series of papers, Witten [11, 12] and Blau and 
Thompson [18] have related 2dYM/BF theory to Chern-Simons theory and G/G WZW models. 
An important upshot of this was in using the Verlinde formula [35] in combination with the 
Riemann-Roch formula to recover the volumes of moduli spaces of flat connections on Riemann 
surfaces: 

dim.M(E fl ) = lim k~* dim *«(E 9 ) dim V g , k . (49) 

The left-hand side was given by the partition function of 2dYM in the topological limit or, 
alternatively, by topological BF theory. Let us recall the definition: 



■^2dYM(S g ,5YM) = 




In the last line we rewrote 2dYM in terms of BF-theory with a quadratic potential. Now the 
point is that in the limit g YM — ► the theory is topological, and the path integral reduces to a 
delta function with support on flat connections. Since we integrate over all gauge connections, 

7 On general Riemann surfaces, there is a second counterterm involving the Euler characteristic. 
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the result is the left-hand side of (49) . Now the right-hand side is given in conformal field theory 
by the modular 5-matrix: 

dimV g , k = Y,S 2 p f 9 • (51) 

A 

It is after all not surprising that this is related to 2dYM if we consider the well-known fact that 
the quantum dimensions reduce to the dimensions of representations of the gauge group in the 
large k limit: 

dim A = lim (52) 

fe^oo opp 

and the partition function of 2dYM on a Ricmann surface of genus g is, in the limit g YM — ► 0: 

^2dYM(S 9 ,0) = ^(dimA) 2 ^ . (53) 

A 

This equation is however only formal, because this limit is in general not well-defined. As we 
already mentioned, in certain cases 2dYM is non-analytic in the coupling constant, and this is 
an indication of the singular nature of the moduli space [18]. 

If we include s punctures with holonomies around them and label them by representations 
Hi, % = 1, . . . , s, the formula for the Verlinde numbers becomes: 

s 

dunV g , s , k = J2sl-p 29 - S U S ^ ■ ( 54 ) 

A i=l 

There is a generalization of the above to finite k due to Blau and Thompson [18]. They 
show that the G/G WZW model on S 9 , which is equivalent to Chern-Simons on S 9 x S 1 , is 
equivalent to BF theory with a compact scalar. Roughly speaking, when we take k — > oo the 
scalar becomes non-compact and we get the topological YM/BF theory back. Obviously, at finite 
k the partition function of this compact BF theory is no longer a delta-function. Nevertheless 
this still computes (51), i.e. the dimensions of the spaces of conformal blocks of the WZW model 
when k — ► oo. 

In the case of the cylinder we have been considering in this paper, the formula (54) is rather 
trivial: it reduces to (S 2 )p jl p 2 = <^ 1/12 , which is simply the satement that the Hilbert space of 
Chern-Simons theory on S x R (the Hilbert space of conformal blocks on S) when £ is a cylinder 
is one-dimensional if [i\ and H2 are conjugate representations, and zero otherwise [28]. 

As we have seen, however, we do not need to consider g YM to be small to get back quantities 
that have an interpretation in the WZW model at finite k. This can be seen as a generalization 
of the relation in [11, 18] between 2dYM and Chern-Simons/WZW to finite k. Indeed, we have 
seen that the partition function of 2dYM computes matrix elements of the type ST q S. Notice 
that the issue of normalization was exactly the same we encountered above for the k — ► oo 
case: 2dYM quantities naturally come normalized as in (52), which are well-defined at large k, 
whereas Chern-Simons naturally computes matrix elements of the type (51), which in particular 
diverge in the classical limit. That is where the normalization (49) comes from. The situation 
for 2dYM with coupling g YM A = ^+g)q ms ^ ea d of zero is exactly the same: we need to normalize 
a matrix element {ST q S)\p/S\ p Spp to be able to compare the theories. 

This seems to suggest that the limit of 2dYM where the imaginary part of the coupling is 
(up to a factor of 2ir) a rational number, and the real part goes to zero, is still some sort of 
topological limit of the theory, precisely analogous to the g YM — ► limit, but not a classical 
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limit. In fact also in this case we can write down a BF theory 

SBF= s4„ Tr * F -&/ E „ d " Tr02 (55> 

corresponding to this YM theory. Notice that, in contrast to [18], the scalar is now non-compact. 

Let us first look at Chern-Simons theory on the sphere, S 3 . This case is rather trivial; as 
mentioned before, it corresponds to 2dYM on the cylinder with trivial representations at the 
two ends and coupling g YM A = In [32] it was found by mirror symmetry that it is described 
by BF theory with a compact scalar and quadratic potential. Notice however that this theory 
lived on a P 1 . This is a natural generalization of the compact BF theory found in [18], which 
lived on E and described Chern-Simons on Ex S 1 , to the case of Chern-Simons on S 3 : one adds 
a quadratic potential, which corresponds to the gluing operator S. In our case, the Yang-Mills 
theory lives on the cylinder, whereas the scalar is non-compact, and the potential is quadratic 
as well. Notice that the gauge coupling is the same as in [32]. It would be interesting to show 
that these theories are the same at the level of the path integral. Notice that our cylinder seems 
to be the same as the cylinder on the B- model side in [32]. On the matrix model side, although 
the measure was unitary, the branes did not go back to themselves after a shift of 2ir around 
the cylinder. 

In the other cases S 3 /Z q , the potential of the BF theory that we obtain is the same but 
divided by 1/q. This is completely analogous to the situation in [32], where the matrix model 
potential had a factor of 1/q. 



5 Discussion and outlook 

Two-dimensional Yang- Mills on a manifold of non-trivial topology is a simple but still interesting 
gauge theory. Besides the interest that it has on its own and its traditional string interpretation 
[2, 6], recent developments mentioned in the introduction [5, 9, 8, 10, 7] indicate that it also 
has useful applications in 4d SYM and topological strings. According to [5], it may even tell us 
something about black holes. 

The topological A-model is effectively described by Chern-Simons theory at large N. If, on 
the other hand, topological amplitudes on certain local threefolds can be described by 2dYM 
[5], one would expect to find direct relations between Chern-Simons theory and 2dYM. It was 
already known that topological 2dYM describes the large k limit of Chern-Simons on manifolds 
of the type E x S 1 and E x R. In this paper we have shown that there is another sector of 2dYM 
on the cylinder - that where the coupling is 2ni times a rational number - which describes 
Chern-Simons on S 3 /Z q at finite k, N. This Yang-Mills theory can be rewritten as a usual BF 
theory with a quadratic potential, and it would be interesting to investigate this theory. 

The S- and T-transformations that we found might also shed some light on the properties 
of the theory at these special values of the coupling. It would be interesting to see whether the 
imaginary part of g Y M c & n be interpreted as a #-angle in YM/BF theory. 

Although we did not take the 't Hooft limit for the coupling, this can still naturally be 
done in our framework as long as k is not much larger than TV, since g YM was proportional to 
7 = FHv- ^his gi yes some hope that the matrix model of 2dYM might have some relevance to 
Chern-Simons in lens spaces. This matrix model will however have different properties from the 
traditional 2dYM matrix model, which generically has a phase transition at finite values of the 
real part of the 't Hooft coupling 8 . In our case, this coupling is purely imaginary. Besides, one 

8 One should however keep in mind the results in [20]. 
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will have to take into account the divergences. Such matrix model might nevertheless still be 
useful, as the natural matrix model that one finds on the Chern-Simons side becomes increasingly 
complicated for larger values of q [32, 33]. Notice that our relation <?ym^ = precisely matches 
the expectation for the topological string, where the effective string coupling is g s = -j^ [32]. It 
would be extremely interesting to see whether the present relation between Chern-Simons and 
2dYM can be obtained geometrically on the topological strings side. We already saw that the 
compact BF theory is obtained by considering branes in the B-model side [32], and it would be 
interesting to show that both theories are equivalent. 

We saw that the Hilbert spaces of this 2dYM and of the WZW model basically agree. The 
periodicity of the exponents with respect to shifts of the coupling by these imaginary values, and 
the fact that we are considering special external states at the two ends of the cylinder, precisely 
conspire so that the infinite sum on the YM side collapses to a sum over the finite Hilbert-space 
of Chern-Simons theory and get the WZW matrix elements. 

An interesting extension of our results would be to include the #-angle in the U(N) case as 
in [5]. This would presumably result in additional linear terms in Frenkel's theorem [14], and 
the generalization should be more or less straightforward. This 0-angle plays a crucial role in 
the relation with topological strings. 

Another interesting extension of this work would be to consider other Riemann surfaces. At 
least naively, we do not really expect the present computations to go through for a Riemann 
surface of arbitrary genus and number of punctures, because adding handles or punctures adds 
factors of dim A, whereas in matrix elements involving products of S and T we expect the 
quantum dimensions. However, we do expect that these ideas can be generalized if we consider 
products of characters. The Boltzmann factors should then determine the value of the Yang-Mills 
coupling. 

In [29], explicit expressions have also been found for Chern-Simons on torus fibrations. The 
expression for the partition function does not involve any external states and it bears some 
resemblance with 2dYM on the torus. It would be interesting to see whether these expressions 
can be related to 2dYM. 

Finally, let us mention that the Brownian motion description of 2dYM made its modular 
transformation properties and some of its analytic structure absolutely explicit and provided us 
with some tools that were useful in relating 2dYM to Chern-Simons theory. In fact we could 
say that Brownian motion densities compute correlators in the WZW model. In [36] it has 
been suggested that Brownian motion could be used to compute certain correlators that involve 
two-matrix models. It would be interesting to explore these connections further. 
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A Lie algebra conventions 

In this appendix we give some of the definitions and conventions used in the paper. 

In this paper we adhere to the nomenclature that calls the affine Weyl chamber a "Weyl 
alcove" , and reserves the name of "chamber" for the finite case. 
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The Weyl alcove C w is defined analogously to the Weyl chamber of the finite group G: 

C w = {\\(w\,ati) > 0,i = 0, ... ,1}, w G W (56) 

where W is the Weyl group, and A is an affine weight. The fundamental Weyl alcove is given 
by taking w = 1. Weights in the fundamental alcove are of the form 

A = A + kuj + rS , (57) 

where k is the level, and with A G Co (the fundamental Weyl chamber), r G R and 5 = (0; 0; 1), 
ljq = (0; 1; 0). P is the weight lattice, P+ the set of dominant (or highest) weights, and -P++ are 
the regular weights, i.e. they are related to those in P + by A = [i + p where A G P++, fJ> G P+. 
P+ are the integrable highest weights, so (A, 6) < k. 

For completeness, we give the explicit expressions for the S- and T-modular matrices used 
in (30). They are represented on affine characters as follows: 

S» = JTT^ \ P /Q V \~ h E eHe"^ (A '"' M) , (58) 
where the central charge is c = k dim (k + g). 

B Divergence of the partition sum 

In this appendix we analyze the leading divergence in (40): 

z(X) = E E e - e ^ va+ ^ V < aV+A >l-T C« va +2« V A) > (59) 

where we already dropped a contribution that vanishes in the limit e — ► 0. The divergence in 
this expression comes from the infinite sum over Q v in the limit e — > 0. In this limit, (59) is 
completely independependent of a v and so the expression is not well-defined. 

Concentrating on the diverging sum over Q v , the relevant e-dependent contribution will be 

f(e,X)= £ e -[^WA] j (60) 

where we shifted A — > A — a v . 
The above can be written as 

f(e, z)=Y, e - e[q2nCn+2qnz] (61) 

where z is linear in the A's and C is the Cartan matrix. This is a generalized theta function 9 
[38]: 



e 



(0\t)q = ^2 ex P [-KT{n + a) n - (n + a) + 2mn ■ b] (62) 



3 We actually use the notation in [37]. 
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It satisfies: 



&[ a b ] (0,r) n = 



r r / 2 VdetU [-a 



(0, l/r) n -i . 



(63) 



So we have 



/(€,«) = e 




1 ~ [-egz/vri 



(eg 2 /^)'V 2 v / detC L 



(0,vr/eg 2 ) c -i , 



(64) 



and we used the fact that the Cartan matrix is non-degenerate. Notice that in the limit e — > 
only the m = term contributes to the function, hence up to finite terms: 



so actually the final expression is independent of z. This is the result we used in (41). 
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